Stability and spectral convergence of Fourier method for nonlinear problems: on the shortcomings of the \(2/3\) de-aliasing method

نویسندگان

  • Claude Bardos
  • Eitan Tadmor
چکیده

The high-order accuracy of Fourier method makes it the method of choice in many large scale simulations. We discuss here the stability of Fourier method for nonlinear evolution problems, focusing on the two prototypical cases of the inviscid Burgers’ equation and the multi-dimensional incompressible Euler equations. The Fourier method for such problems with quadratic nonlinearities comes in two main flavors. One is the spectral Fourier method. The other is the 2/3 pseudo-spectral Fourier method, where one removes the highest 1/3 portion of the spectrum; this is often the method of choice to maintain the balance of quadratic energy and avoid aliasing errors. Two main themes are discussed in this paper. First, we prove that as long as the underlying exact solution has a minimal C1+α spatial regularity, then both the spectral and the 2/3 pseudo-spectral Fourier methods are stable. Consequently, we prove their spectral convergence for smooth solutions of the inviscid Burgers equation and the incompressible Euler equations. On the other hand, we prove that after a critical time at which the underlying solution lacks sufficient smoothness, then both the spectral and the 2/3 pseudo-spectral Fourier methods exhibit nonlinear instabilities which are realized through spurious oscillations. In particular, after shock formation in inviscid Burgers’ equation, the total variation of bounded (pseudo-) spectral Fourier solutions must increase with the number of increasing modes and we stipulate the analogous situation occurs with the 3D incompressible Euler equations: the limiting C. Bardos University of Paris 7, Denis Diderot, Laboratory Jacques Louis Lions, University of Paris 6, Paris, France e-mail: [email protected] E. Tadmor (B) Department of Mathematics, Center of Scientific Computation and Mathematical Modeling (CSCAMM), Institute for Physical sciences and Technology (IPST), University of Maryland, College Park, MD 20742-4015, USA e-mail: [email protected]

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Error Estimation of Pseudospectral Method for Solving the Barotropic Vorticity Equation

The barotropic vorticity equation model is very important in the research of meteorological science and applied mathematics. Many scientists pay attention to the research of numerical methods of this equation. The early works were mainly concerned with finite difference methods [2, 4, 8, 9, 13]. The spectral method has a convergence rate of ”infinite” order, i.e., the error decays faster than a...

متن کامل

Stability and Convergence Analysis of Fully Discrete Fourier Collocation Spectral Method for 3-D Viscous Burgers' Equation

This paper analyzes the stability and convergence of the Fourier pseudospectral method coupled with a variety of specially designed time-stepping methods of up to fourth order, for the numerical solution of a three dimensional viscous Burgers’ equation. There are three main features to this work. The first is a lemma which provides for an L2 and H 1 bound on a nonlinear term of polynomial type,...

متن کامل

A numerical technique for solving a class of 2D variational problems using Legendre spectral method

An effective numerical method based on Legendre polynomials is proposed for the solution of a class of variational problems with suitable boundary conditions. The Ritz spectral method is used for finding the approximate solution of the problem. By utilizing the Ritz method, the given nonlinear variational problem reduces to the problem of solving a system of algebraic equations. The advantage o...

متن کامل

Solving a nonlinear inverse system of Burgers equations

By applying finite difference formula to time discretization and the cubic B-splines for spatial variable, a numerical method for solving the inverse system of Burgers equations is presented. Also, the convergence analysis and stability for this problem are investigated and the order of convergence is obtained. By using two test problems, the accuracy of presented method is verified. Additional...

متن کامل

Fourth-order numerical solution of a fractional PDE with the nonlinear source term in the electroanalytical chemistry

The aim of this paper is to study the high order difference scheme for the solution of a fractional partial differential equation (PDE) in the electroanalytical chemistry. The space fractional derivative is described in the Riemann-Liouville sense. In the proposed scheme we discretize the space derivative with a fourth-order compact scheme and use the Grunwald- Letnikov discretization of the Ri...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:
  • Numerische Mathematik

دوره 129  شماره 

صفحات  -

تاریخ انتشار 2015